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Closed, singularity- free, inflationary cosmological models have recently been studied in the context 
of general relativity. Despite their appeal, these so called emergent models suffer from a number of 
limitations. These include the fact that they rely on an initial Einstein static state to describe the 
past eternal phase of the universe. Given the instability of such a state within the context of general 
relativity, this amounts to a very severe fine tuning. Also in order to be able to study the dynamics 
of the universe within the context of general relativity, they set the initial conditions for the universe 
in the classical phase. Here we study the existence and stability of such models in the context of 
Loop Quantum Cosmology and show that both these limitations can be partially remedied, once 
semi-classical effects are taken into account. An important consequence of these effects is to give 
rise to a static solution (not present in GR), which dynamically is a centre equilibrium point and 
located in the more natural semi-classical regime. This allows the construction of emergent models 
in which the universe oscillates indefinitely about such an initial static state. We construct an 
explicit emergent model of this type, in which a non-singular past eternal oscillating universe enters 
a phase where the symmetry of the oscillations is broken, leading to an emergent inflationary epoch, 
while satisfying all observational and semi-classical constraints. We also discuss emergent models in 
which the universe possesses both early- and late-time accelerating phases. 
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I. INTRODUCTION 

One of the fundamental questions of modern cosmol- 
ogy is whether the universe had a definite origin or 
whether it is past eternal. The central paradigm for 
structure formation in the universe is the inflationary 
scenario. (For a review, see, e.g., Ref. pj). Under very 
general conditions, inflation is future eternal, in the sense 
that once inflation has started, most of the volume of the 
universe will remain in an inflating state |2|. Given cer- 
tain assumptions, however, it has been argued Q that 
inflation can not be past null complete - and therefore 
past eternal - if it is future eternal. 

Leading alternatives to inflation that are motivated by 
recent advances in string/M-theory are the pre-big bang 
0,01 and ekpyrotic/cyclic £ 6j scenarios, respectively. The 
fundamental postulate of the pre-big bang model is that 
the set of initial data for the universe lies in the infinite 
past in the perturbative regime of small string coupling 
and spacetime curvature The universe then evolves 
into a strongly coupled and highly curved regime before 
exiting into the standard, post-big bang phase. 

On the other hand, the big bang singularity is inter- 
preted in the cyclic scenario in terms of the collisions of 
two co-dimension one branes propagating in a higher- 
dimensional spacetime 0. In this picture, the branes 
undergo an infinite sequence of oscillations where they 
move towards and subsequently away from each other. 
From the perspective of a four-dimensional observer, the 
collision of the branes is interpreted as the bounce of a 
collapsing universe into a decelerating, post-big bang ex- 
pansion. 

Despite these attractive properties, however, the pro- 
cess that leads to a non-singular transition between the 
pre- and post-big bang phases is unclear. The search for 



singularity-free inflationary models within the context of 
classical General Relativity (GR) has recently led to the 
development of the emergent universe scenario 

EES In 

this model, the universe is positively curved and initially 
in a past eternal Einstein static (ES) state that eventu- 
ally evolves into a subsequent inflationary phase. Such 
models are appealing since they provide specific examples 
of non-singular (geodesically complete) inflationary uni- 
verses. Furthermore, it has been proposed that entropy 
considerations favour the ES state as the initial state for 
our universe [jj. 

Classical emergent models, however, suffer from a 
number of important shortcomings. In particular, the 
instability of the ES state (represented by a saddle equi- 
librium point in the phase space of the system) makes 
it extremely difficult to maintain such a state for an in- 
finitely long time in the presence of fluctuations, such as 
quantum fluctuations, that will inevitably arise. More- 
over, the initial ES state must be set in the classical do- 
main in order to study the dynamics within the context 
of GR. A more natural choice for the initial state of the 
universe would, however, be in the semi-classical or quan- 
tum gravity regimes. 

A leading framework for a non-perturbative theory of 
quantum gravity is loop quantum gravity (LQG). (For re- 
views, see, e.g., Ref. |lfi ITll|L Recently, there has been 
considerable interest in loop quantum cosmology (LQC), 
which is the application of LQG to symmetric states |12| . 
Within the framework of LQC, there can exist a 'semi- 
classical' regime, where spacetime is approximated by a 
continuous manifold, but where non-perturbative quan- 
tum corrections modify the form of the classical Einstein 
field equations |l3j|. 

Motivated by the above discussion, the present paper 
studies the existence and nature of static solutions within 
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the framework of semi-classical LQC in the presence of a 
self-interacting scalar field. Significant progress has re- 
cently been made in understanding the dynamics of this 
system, and its relevance to inflation and non-singular 
behaviour [3, EH [l^, [T3, EH E3, HH - In addition 
to the ES universe of classical cosmology, LQC effects 
result in a second equilibrium point in the phase space 
that we refer to as the 'loop static' (LS) solution. Cru- 
cially, this solution corresponds to a centre equilibrium 
point in the phase space and, consequently, the universe 
may undergo a series of (possibly) infinite, non singular 
oscillations about this point. During these oscillations, 
the scalar field can be driven up its potential ^3- This 
leads us to consider a new picture for the origin of the 
universe, where the universe is initially oscillating about 
the LS static solution in the infinite past and eventu- 
ally emerges into a classical inflationary era. The model 
shares some of the attractive features of the pre-big bang 
and cyclic scenarios, in the sense that it is past eternal, 
although it exhibits a significant difference in that the cy- 
cles are broken by inflationary expansion. Additionally, 
unlike the pre-big bang and cyclic scenarios, the emer- 
gent model considered here is genuinely non-singular. 

The paper is organised as follows. The existence and 
nature of static universes in semi-classical LQC is stud- 
ied in Section 2. Section 3 discusses the dynamics that 
leads to the emergence of an inflationary universe and 
a particular inflationary model that is consistent with 
present-day cosmological observations is outlined in Sec- 
tion 4. In Section 5 we discuss a class of models where 
the field that drives the pre-inflationary oscillations may 
also act as the source of dark energy in the present-day 
universe. We conclude with a discussion in Section 6. 



II. STATIC UNIVERSES 



field. Equations and J5J together admit the first in- 
tegral (the Friedmann equation): 
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where as usual K = 0,±1 parametrises the spatial cur- 
vature. Combining Eq. JSJ) with Eq. JJJ yields 
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The Einstein static universe is characterised by the 
conditions K = 1 and a = = a. In the presence of 
a scalar field with a constant potential V, it is straight- 
forward to verify that there exists a unique solution sat- 
isfying these conditions, where the scale factor is given 
by a = ao with ao = (iirlp^ 2 ) 1 / 2 . There are two im- 
portant points to note regarding this solution. Firstly, it 
represents a saddle equilibrium point in the phase space, 
which is unstable to linear perturbations. Secondly, in 
order to enable the analysis to be performed within the 
context of GR, it is necessary to assume that ao lies in 
the classical domain. 

In the emergent universe scenario, it is assumed that 
the potential becomes asymptotically flat in the limit 
(j) — > — oo and that the initial conditions are specified 
such that the ES configuration represents the past eter- 
nal state of the universe, out of which the universe slowly 
evolves into an inflationary phase. However, the insta- 
bility of the ES universe ensures that any perturbations 
- no matter how small - rapidly force the universe away 
from the ES state, thereby aborting the scenario. 

In the next subsection we shall see that employing a 
more general LQC setting can partially resolve both these 
issues. 



A. Einstein static universe in classical gravity 



B. Static solutions in semi-classical LQC 



Before considering static solutions within the context 
of LQC, it is instructive to review the corresponding re- 
sults for classical Einstein gravity. 

Throughout this paper, we consider an isotropic and 
homogeneous Friedmann-Robertson- Walker (FRW) uni- 
verse sourced by a scalar field with energy density and 
pressure given by p = \ (j) 2 + V{<j>) and p = | <fr 2 — V (</>), 
respectively, where V(<f>) represents the self-interaction 
potential of the field. The Raychaudhuri equation for 
such a universe is given by 
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and local conservation of energy-momentum implies the 
scalar field satisfies 

4> + 3H<j) + V'(<t>) = , (2) 

where a prime denotes differentiation with respect to the 



It has recently been shown that in the semi-classical 
regime the cosmological evolution equations become 
modified When restricted to FRW backgrounds, this 
regime is characterised by the scale factor of the universe 
lying in the range a, < a < a*, where a; = y^yZpi deter- 
mines the scale for the discrete quantum nature of space- 
time to become important and 7 ~ 0.274 is the Barbero- 
Immirzi parameter |22j. and a» = \J 7i/3Zpi, where j is 
a parameter that arises due to ambiguities in the quanti- 
zation procedure j^. It must ta ke p ositive, half-integer 
values but is otherwise arbitrary 38] . The standard clas- 
sical cosmology is recovered above the scale a* and the 
parameter j therefore sets the effective quantum grav- 
ity scale. The dynamics in this regime is determined 
by replacing the inverse volume term a -3 that arises in 
the classical matter Hamiltonian = ^a~ 3 p 2 ^ + a 3 V((j)) 
with a continuous function dj(a) that approximates the 
eigenvalues of the geometrical density operator in LQC. 
This 'quantum correction' function is given by dj(a) = 
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D(q)a 3 , where 

( 8 \ 3 r_ r / , ii ,ii 

D(q) = f-J g a|7[(g + l)T _|g_l|T 



(q+l)i - sgn(q - l)\q - 1\ 7 -]Y , (5) 



and q = (a/a*) 2 . As the universe evolves through the 
semi-classical phase, this function varies as D a a 15 for 
a <C a*, has a global maximum at a ss a*, and falls 
monotonically to D = 1 for a > a*. 

The effective field equations then follow from the 
Hamiltonian fl3l: 



H = - 
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where p^, — d 3 ■ 1 cf) is the momentum canonically conjugate 
to the scalar field. The Raycharduri equation becomes 
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3L> V 4 dma / 3 
and the modified scalar field equation takes the form 
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The first integral of Eqs. JIJ) and JHJ is given by the 
modified Friedmann equation 
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where in the LQC context K can only take values 0, +1. 
Combining Eq. Q with Eq. 10 then implies that 
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Equation l(7|l- (|10|l can be rewritten in the form of a 
three-dimensional dynamical system in terms of variables 
{a, if, V } H3 ■ (The present study extends the results of 
[2fJ and, more importantly, provides an analytical ac- 
count of the dynamics). Employing the Friedmann equa- 
tion © to eliminate the scalar field's kinetic term, and 
assuming that dV/dip can be expressed as a function of 
the potential, allows the complete dynamical system to 
be described by equations: 



H = (Sirl^V - 3H 2 ) 1 - - 
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where the function A(a) = din D /din a. In this Section, 
with the aim of considering emergent universes we shall 
take K = +1 and initially consider constant potentials. 
In this case, Eq. (|13fl is trivially satisfied and the sys- 
tem reduces to the two-dimensional autonomous system 

The static solutions (d = = d) then correspond to 
the equilibrium points of this system, which are given by 



H cq — 0, 
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where et cq is given by solutions to the constraint equation 



A(a cq ) = B(a oq ), B(a) 
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It follows from Eq. (|1(J[) that the field's kinetic energy at 
the equilibrium points is given by 
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Equation (|15fl implies that a necessary and sufficient 
condition for the existence of the equilibrium points is 
that the functions A(a) and B(a) should intersect. 

Here we are interested in how the properties of the 
equilibrium points alter as the potential is varied. This 
is best seen by considering how the functions A(a) and 
B(a) (and hence their points of intersection) change. The 
function A has the form shown in the top panel of Fig. ^ 
It asymptotes to the constant value A = 15 at a = 0, 
decreases to a minimum value of A m i n = —5/2 at a = 
a», and then asymptotes to zero from below as a — * oo. 
Increasing the parameter j increases the value of a* and 
this results in moving the turning point of the function 
A to larger values of the scale factor. The important 
point to note, however, is that the qualitative form of 
the function A remains unaltered for all values of j, as 
can be seen from Fig. ^ 

The qualitative nature of the function B, on the other 
hand, is sensitive to the sign of the potential, as can be 
seen from the bottom panel of Fig. We briefly discuss 
the behaviour of B for each case in turn. For V = 0, 
it is given by the (solid) horizontal line B = 4. For 
V > it represents a hyperbola (dot-dashed curve) with 
a single (since the scale factor a is non-negative) vertical 
asymptote given by 
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The region to the left of this asymptote defines the re- 
gion in the {A/B, a} plane where the reality condition, 
cf> 2 > 0, is satisfied. Thus, the upper-right branch of the 
hyperbola plays no role in determining the existence of 
the equilibrium points. In the limit of a — > 0, the func- 
tion B — > 4 which coincides with the value of the function 
in the case V = 0. As a is increased, the function takes 
progressively smaller values as the vertical asymptote is 
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approached. As V tends to zero from above, it causes 
the asymptote to move to progressively larger values of 
a, and we may therefore formally view the V = case as 
the limit where the asymptote moves to infinity. Finally, 
for negative potentials the qualitative behaviour of the 
function B is changed to the dashed line in the bottom 
panel of Fig. and there is no vertical asymptote. As in 
the case of positive potentials, B — > 4 as a — > 0, but the 
function B now increases monotonically as a increases, 
ultimately tending to the value 6 as a — > oo. As the 
potential tends to zero from below, the function B still 
tends to the asymptotic value B = 6, but at larger a. 
An important point to note is that for all values of V, 
the function B satisfies the condition B < 6 (for physi- 
cally relevant regions of the {A/B,a} plane), and hence 
equilibrium points can only occur when A = B < 6. 

Once the positions of the equilibrium points have been 
determined, their nature can be found by linearising 
about these points. The eigenvalues are given by 
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When A 2 < 0, this leads to imaginary eigenvalues and 
implies that the equilibrium point is a centre, whereas 
the point is a saddle when A 2 > 0. Although this ex- 
pression is rather complicated, it can be verified that for 
A < 6 (which is a necessary condition for the existence of 
equilibrium points), A 2 is negative for a < a* and posi- 
tive for a > a*. Hence we have the important result that, 
equilibrium points occurring in the semi- classical regime 
are centres and those occurring in the classical regime are 
saddles. 

To determine the existence of equilibrium points, we 
again consider the different signs of the potential sepa- 
rately. 

V = 0: In the case of a massless scalar field with V = 0, 
the condition (|15|l has a single solution given by A = 4, 
implying a single equilibrium point. The position of this 
equilibrium point can be seen from the intersection point 
in the top panel of Fig. [21 Since this point occurs at 
a < a», it is a centre, with its phase portrait consisting 
of a continuum of concentric orbits (see the bottom panel 
of Fig. |3. 

V > 0: As described above, for small positive values 
of the potential, the vertical asymptote of the function 
B is far to the right of the origin, which results in two 
points of intersections between the functions A and B 
(see the top panel of Fig. [3J. There are therefore two 
equilibrium points in this case: the first occurs at a < a* 
and hence is a centre, while the second occurs at a > a* 
and is therefore a saddle. As the potential is increased, 
the asymptote moves towards the origin causing the equi- 
librium points to eventually coalesce when V = V*. This 
occurs at the point of tangency of the functions A and 
B, i.e. at the minimum (which is a cusp) of the function 
A located at a = a*. Thus, using B and the fact that 
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FIG. 1: The top panel shows the plots for the function A 
against a for j = 100 (solid line), j = 1000 (dashed line) 
and j = 3000 (dot-dashed line). The bottom panel shows 
plots of the function B against a for V — (solid line) , V > 
(dot-dashed line) and V < (dashed line) . Note the dramatic 
change in the shape of the function as the sign of V is altered. 
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For values of V > V*, the asymptote moves further 
towards the origin. Consequently, the functions A and 
B will no longer intersect in the part of the plane which 
satisfies the reality condition, and therefore there will be 
no equilibrium points in this case. It is also easy to un- 
derstand the effect of changing the parameter j on V*. 
Increasing j increases a* , which means that the point of 
tangency occurs for smaller values of the potential. This 
implies that larger values of the quantization parameter, 
and hence a*, lead to smaller values of the potential be- 
yond which no equilibrium points can exist. 
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FIG. 2: The top panel shows the plots of functions A (solid 
line) and B (dashed line) against a for V = 0. The vertical 
dotted line denotes the position of a = a*. The bottom panel 
is the corresponding phase portrait demonstrating the centre 
equilibrium point that occurs in this case. The plot is com- 
pactified using x = arctan(iif) and y — arctan(lna), in order 
to present the entire phase space. The vertical dotted lines 
in the bottom panel demarcate the region in which the Hub- 
ble parameter is less than the Planck scale. The horizontal 
dotted line marks a», where the quantum regime begins. All 
axes are in Planck units, and j — 100. 



FIG. 3: The top panel shows the plots of functions A (solid 
line) and B (dashed line) against a for V = 0.005. For this 
value of V, the equilibrium points are close to coalescing. The 
vertical dotted line denotes the position of a — a«. The bot- 
tom panel is the corresponding phase portrait demonstrating 
the centre equilibrium point and saddle point which occurs in 
this case with the dashed line indicating the separatrix. The 
shaded area is the unphysical region where if) 2 would be neg- 
ative. The plot is compactified as in Fig. H The dotted lines 
in the bottom panel are as described in Fig. All axes are 
in Planck units, and j = 100. 



V < 0: Negative potentials are known to arise in 
string/M-theory compactifications (e.g., |24|) and are of 
interest in connection with the ekpyriotic/cyclic mod- 
els considered recently |6(. In such cases the qualitative 
change in the form of B results in only one equilibrium 
point, which occurs inside the semi-classical region and 
hence is a centre (see the bottom panel of Fig. 0J. Fur- 
thermore, the centre equilibrium point will persist for any 
negative value of V. 

To summarise, the important consequence of consider- 
ing LQC effects in this context is that they admit two 
possible static solutions, rather than the single solution 



in the case of GR. The first static solution corresponds 
to a saddle point (as in GR) and is referred to as the Ein- 
stein Static (ES) solution. The second static solution is a 
direct consequence of LQC effects. It is a centre and we 
shall refer to it as the loop static (LS) solution. We have 
shown analytically that the LS solution arises for a wide 
range of values for the potential given by V € (— Vw, V*), 
where the lower bound —Vn is determined by the need 
to satisfy the Planck bounds. The importance of the LS 
solution is that, in contrast to the ES solution present in 
GR, slight perturbations do not result in an exponential 
divergence from the static universe, but lead instead to 
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FIG. 4: The top panel shows the plots of functions A (solid 
line) and B (dashed line) against a for V = —0.013. The ver- 
tical dotted line denotes the position of a — a*. The bottom 
panel is the corresponding phase portrait demonstrating the 
centre equilibrium point that occurs in this case. The plot 
is compactified as in Fig. [5] The dotted lines in the bottom 
panel are as described in Fig. [5] All axes are in Planck units 
and j = 100. 



oscillations about it. Moreover, since the LS solution al- 
ways lies within the semi-classical region a < a* , it is ide- 
ally suited to act as the past asymptotic initial state for 
an emerging universe. Indeed, it is the only equilibrium 
point in models with a vanishing or negative potential. 



In the following Section, we employ these novel features 
of LQC dynamics within the context of the emergent in- 
flationary universe. 



III. EMERGENT INFLATIONARY UNIVERSE 
IN LQC 

A. The Dynamics of Emergence 

The phase plane analysis of the previous Section de- 
termined the qualitative LQC dynamics for the case of 
a constant potential. However, any realistic inflationary 
model clearly requires the potential to vary as the scalar 
field evolves. Nevertheless, the constant potential dy- 
namics (see Figs. provides a good approximation to 
more general dynamics if variations in the potential are 
negligible over a few oscillations. This implies that cos- 
mic dynamics with a variable potential may be studied by 
treating the potential as a sequence of separate, constant 
potentials. Here, with the emergent inflationary models 
in mind (see below), we shall consider a general class of 
potentials that asymptote to a constant V_oc < V* as 
(j) — > — oo and rise monotonically once the value of the 
scalar field exceeds a certain value. 

As we saw in the previous Section, for any constant po- 
tential V <V*, there exists a region of parameter space 
where the universe undergoes non-singular oscillations 
about the point LS. For V > V*, the equilibrium points 
LS and ES merge and the phase plane then resembles that 
of classical GR: a collapsing universe bounces and asymp- 
totically evolves into a de Sitter (exponential) expansion 
in the infinite future. This leads us to propose the follow- 
ing picture for the origin of the universe. The universe is 
initially at, or in the neighbourhood of, the static point 
LS, with the field located on the plateau region of the 
potential with a positive kinetic energy, <p- m it > 0. The 
universe undergoes a series of non-singular oscillations 
in a (possibly) past-eternal phase with the field evolv- 
ing monotonically along the potential. As the magnitude 
of the potential increases, the cycles are eventually bro- 
ken by the emergence of the universe into an inflationary 
epoch. 

We proceed to discus the dynamics of this scenario in 
more detail. The universe will exhibit cyclic behaviour 
around the LS point for a very wide range of initial val- 
ues </>i„it when V-oo < 0. If < V-oo < V*, the range 
becomes more limited as V-oo is increased. An impor- 
tant feature to note is that in all cases the field's kinetic 
energy never vanishes during a given cycle. In the case 
of a positive potential, the unphysical region of phase 
space where the field violates the null energy condition 
lies outside the cyclic region that is bounded by the sep- 
aratrix. Indeed, for a constant potential, the evolution 
of the scalar field is determined by integrating Eq. ©: 



"hut 



D 



(20) 



Thus, the field will vary monotonically along the poten- 
tial and eventually reach the region where the potential 
begins to rise. 

Increasing the magnitude of the potential over a series 
of cycles has the effect of moving the location of the LS 
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point to progressively higher values of the scale factor, 
although the shift is moderate. In the case of a posi- 
tive potential, for example, a necessary condition for the 
existence of LS is that 4 < A(a cq ) < —2.5, where the 



function A is defined in Eq. (|15|l . and this corresponds 
to the range 0.94a* < a < a*. Eq. (|16|l then implies 
that the field's kinetic energy does not alter significantly, 
since the universe remains in the vicinity of LS. On the 
other hand, the saddle point ES occurs at progressively 
smaller values of the scale factor as the magnitude of 
the potential increases. As discussed in Section 2.2, this 
equilibrium point occurs in the range a* < a < oo, where 
the limits are approached as V — > V* and V — ► 0, respec- 
tively. 

The overall effect of increasing the potential, therefore, 
is to distort the separatrix in the phase plane, making it 
narrower in the vertical direction but introducing little 
change to the position of the LS point. This implies that 
the dynamics varies only slightly from cycle to cycle for 
orbits that are close to the LS point. If the magnitude 
of the potential continues to increase as the field evolves, 
however, a cycle is eventually reached where the trajec- 
tory that represents the universe's evolution now lies out- 
side the finite region bounded by the separatrix and this 
effectively breaks the oscillatory cycles. From a physi- 
cal point of view, the magnitude of the field's potential 
energy, relative to its kinetic energy, is now sufficiently 
large that a recollapse of the universe is prevented, i.e., 
the strong energy condition of GR is violated, thereby 
leading to accelerated expansion. The field decelerates as 
it moves further up the potential, subsequently reaching 
a point of maximal displacement and then rolling back 
down. If the potential has a suitable form in this region, 
slow-roll inflation will occur. 



B. The Energy Scale of Inflation 

An important question to address is the energy scale at 
the onset of inflation, Vinf . In general, inflation begins (in 
the classical regime) when the strong energy condition is 
violated: 



V(</>inf) ~ 4>t 



,f ) 



(21) 



and, moreover, the structure of the phase space indicates 
that the potential energy of the field remains dynami- 
cally negligible until the onset of inflation [3jJ. In this 
case, the evolution of the scalar field prior to inflation is 
determined by Eq. l(2U|) . 

The inflationary energy scale may then be estimated 
by considering the penultimate cycle before the onset of 
inflation. The turnaround in the expansion occurs when 
the field's energy density satisfies p = 3/(87d P1 a 2 ). Since 
the energy density of the field will not have changed sig- 
nificantly at the equivalent stage of the following cycle, 
the scale factor at the onset of inflation is determined 
approximately by a- ln { ~ (47rZ P1 Vhif) -5 . Substituting this 



condition into Eq. I|20[l then yields an estimate for the 
magnitude of the potential in terms of initial conditions: 
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For a universe located near to the equilibrium point 
LS, the scale factor is given by ai n it = /&*, where 0.94 < 
/ < 1 and it may be further verified that Mnit =0(1) 
in this range. The field's initial kinetic energy is then 
determined by Eq. JTSJ): <j)f nit « 3/(47r/|, 1 a 2 q ). It follows 
from Eq. I|22|) . therefore, that a universe 'emerging' from 
the semi-classical LQC phase near to the LS point will 
begin to inflate when 
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As expected, this is in good agreement with the neces- 
sary condition 1)19(1 for the coalescence of the equilibrium 
points LS and ES, since the scale factor can not evolve 
until V > V* if the universe is located on or near the LS 
point. 

A precise measure on the set of initial data is presently 
unknown, and we should therefore consider other possible 
initial conditions. Having investigated the regime ai n i t ~ 
a* , a second possibility is to consider initial conditions 
where ~ aj n ;t <C a*. In this regime, the quantum 
correction factor JSJ asymptotes to a power-law, D w 
(12/7) 6 (a/a») 15 , and Eq. (22Jl is then equivalent to the 
condition 



20 



12 



Vinf 

4 

m pl 



'Pi • 



(24) 



where we have defined the ratio /3 = a; n it/&*. This ratio 
may be constrained by imposing two necessary conditions 
for the semi-classical framework to be valid. Firstly, the 
initial conditions should be set in the regime where space- 
time is approximated by a smooth manifold, ai/a; n it < 1, 
and this leads to the lower limit: 



0> 



(25) 



Secondly, the scalar field's kinetic energy must not exceed 
the Planck scale at the onset of the classical regime. Since 
the anti-frictional effects in the modified field equation 
(jHJ) accelerate the field when a < a*, the field's kinetic 
energy must initially be sub-Planckian, |0i n it|/?7ipi < 1- 
Eq. Q24H then leads to an upper limit on p: 



(3< 



Vinf 
4 

m P1 



(26) 



and combining the constraints l|25|) and (|26|l results in a 
necessary condition for the onset of inflation when aj n i t <C 



m P1 



2 



(27) 
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FIG. 5: The figure depicts two possible forms of emergent 
potentials that allow for conventional re-heating. The solid 
line is the form of the potential motivated by the inclusion of 
a R 2 term in the Einstein-Hilbert action, and the dotted line 
the form motivated by the inclusion of higher-order terms. 



In other words, for a given value of the parameter j, 
inflation must begin above the scale (2/j)"2"mp 1 . 

Conditions l|23|) and (|27fl both imply that the infla- 
tionary energy scale is higher for lower values of the pa- 
rameter j. Indeed, it is comparable to the Planck scale 
for j < O(10) and this generic behaviour is only weakly 
dependent on the initial conditions. This has implica- 
tions for the asymptotic form of the potential as the field 
reaches progressively higher values. Unless the param- 
eter j is sufficiently large, it is unlikely that the oscil- 
latory dynamics will end if the potential asymptotes to 
a constant value, or reaches a local maximum, that is 
significantly below the Planck scale. In this sense, there- 
fore, the scenario outlined above favours potentials that 
increase monotonically once the value of the scalar field 
has exceeded some critical value. 

In the following Section, we consider a concrete exam- 
ple of a potential that exhibits the appropriate asymp- 
totic behaviour. 



IV. A SPECIFIC MODEL OF AN EMERGING 
UNIVERSE 

From a dynamical point of view, the emergent uni- 
verse scenario can be realised within the context of semi- 
classical LQC if the potential satisfies a number of rather 
weak constraints. Asymptotically, it should have a hor- 
izontal branch as <j> — > — oo such that dV/d<j> — > and 
increase monotonically in the region <f> > grO w, where 
without loss of generality we may choose </> gr0 w = 0. The 
reheating of the universe imposes a further constraint 
that there should be a global minimum in the potential 
at Knin = if reheating is to proceed through coherent 



oscillations of the inflaton. Since inflation will end after 
the field has rolled back down the potential, this should 
occur for <j> < grow . Finally, the region of the potential 
that drove the last 60 e-foldings of inflationary expan- 
sion is then constrained by cosmological observations, as 
in the standard scenario. 

Examples of potentials that exhibit these generic prop- 
erties are illustrated in Fig. |3J It is interesting that 
potentials of this form have been considered previously 
in a number of different settings, including cases where 
higher-order curvature invariants of the form 



L 



1 



JV 



N 



rl l — l 



6i R l 



(28) 



are introduced into the Einstein-Hilbert action, where R 
is the Ricci scalar, ej are coupling constants and t\ = 1. 
Such corrections are required when attempting to renor- 
malize theories of quantum gravity pEj and also arise 
in low-energy limits of superstring theories |2(j| . In gen- 
eral, such theories are conformally equivalent to Einstein 
gravity sourced by a minimally coupled, self-interacting 
scalar field. For example, potentials with a nonzero 
asymptote at — > — oo (as shown by the solid line in 
Fig. [3J can be obtained from theories that include a 
quadratic term in the action, whereas those with a zero 
asymptotic value and a local maximum (shown by a dot- 
ted line) arise when cubic and higher-order terms in the 
Ricci scalar are introduced [53, • In general, all these 
potentials possess a global minimum at V — 0. Finally, 
potentials of this form have also been considered within 
the context of the classical emergent universe Q . 

Motivated by the above discussion, we consider, as an 
example, the potential 



V = a Ueiq>(/3<f>/V3) - 1 



(29) 



where a and (3 are constants. This potential is qual- 
itatively similar to that illustrated by the solid line in 
Fig. [S] The parameters in the potential are constrained 
primarily by the Cosmic Microwave Background (CMB) 
anisotropy power spectrum. Assuming inflation proceeds 
in the region <p 3> 0, the parameter f3 determines the 
spectral index of the scalar perturbation spectrum to- 
gether with the relative amplitude of the gravitational 
wave perturbations. The parameter a is then constrained 
by the COBE normalization of the power spectrum on 
large-scales [29]. 

We chose a = 10~ 12 and (3 = 0.1 as representing typi- 
cal values satisfying the constraints imposed recently by 
the WMAP satellite |30| and numerically integrated the 
field equations 0-© for a universe starting from an ini- 
tial state close to the LS static state. The results of the 
integration are illustrated in Figs. HO and [7| The field 
starts in the asymptotically flat region of the potential 
and gradually increases in value, as the scale factor os- 
cillates about the LS point. The field moves past the 
minimum and climbs up the potential. The scale fac- 
tor continues to oscillate until the field reaches the point 
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FIG. 6: Time evolution of the scale factor (top panel) and 
scalar field (bottom panel) with the field initially on the 
asymptotically flat region of the potential I29H with a — 10~ 12 
and f3 — 0.1. The field increases in value from initial condi- 
tions close to the LS (centre) solution defined by Eq. J14I . 
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FIG. 7: Plots illustrating the magnification of Fig. |S|around 
the region where emergence commences, the oscillations cease 
and inflation begins. 



where it slows down significantly, thereby bringing the 
oscillations to an end and initiating the inflationary ex- 
pansion. The behaviour at this stage is qualitatively sim- 
ilar to that discussed previously for a quadratic potential 

m 



V. EMERGING QUINTESSENTIAL INFLATION 

One drawback of reheating through infiaton decay in 
the emergent scenario is that the coupling of the scalar 
field to the standard model degrees of freedom must be 



strongly suppressed if the field is to survive for a (possi- 
bly) infinite time as it emerges from the oscillatory semi- 
classical phase. It is more natural, therefore, to invoke a 
'sterile' infiaton that is not coupled directly to standard 
model fields, and where reheating proceeds through an 
alternative mechanism such as gravitational particle pro- 
duction In this case, the potential need not exhibit 
a minimum and the field could continue to roll back down 
the potential towards <f> — > —oo at late times. 

It is notable that the general requirements for a ster- 
ile infiaton with a potential exhibiting a decaying tail as 
<j> — > — oo are precisely those features that are character- 
istic of the quintessential inflationary scenario, where the 
field that drove early universe inflation is also identified 
as the source of dark energy today j^. In the present 
context, this suggests that the scalar field could play a 
three-fold role in the history of the universe, acting as 
the mechanism that enables the universe to emerge into 
the classical domain, and then subsequently driving both 
the early- and late-time accelerated expansion. 

The specific constraints that must be satisfied by the 
potential in standard quintessential inflation have been 
considered in detail in Ref. |33j . In particular, one 
of the simplest asymptotic forms for the low-energy 
tail that simultaneously leads to tracking behaviour and 
satisfies primordial nucleosynthesis bounds is given by 
V oc (m/cf>) k exp(A0/mpi), where m, k and A are con- 
stants. Cyclic behaviour in LQC will arise for such a 
potential. 

A further requirement is that the potential must be 
sufficiently steep immediately after the end of inflation if 
the field's energy density is to redshift more rapidly than 
the sub-dominant radiation component. This requires a 
second point of inflexion in the potential, as illustrated 
qualitatively in Fig. [SJ Beyond this region, the poten- 
tial must continue to rise in order for the oscillatory dy- 
namics to come to end and, as discussed above, this is 
expected to occur near to the Planck scale. From a dy- 
namical point of view, there are no further constraints 
on how rapidly the potential energy need increase in this 
region. The only remaining consideration is that a phase 
of successful slow-roll inflation should arise as the field 
rolls back down the potential. Given the ease with which 
the infiaton is able to move up the potential due to LQC 
effects, we anticipate that any potential satisfying the 
existing constraints for successful quintessential inflation 
will also lead to a successful emergence of the classical 
universe. 



VI. DISCUSSION 

In this paper, we have investigated the occurrence 
of static solutions in loop quantum cosmology settings 
sourced by a scalar field with a constant potential. We 
have shown that there are in principle two such solutions, 
depending upon the sign and magnitude of the potential. 
The point ES is always in the classical domain and cor- 
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FIG. 8: Plots illustrating the generic form of the potential 
that leads to early- and late-time accelerating phases. The 
potential exhibits a decaying tail as <j> — > — oo. As the field 
moves up this tail and increases in value, the universe can 
oscillate about the LS point. The inflationary regime rises 
(possibly) towards the Planck scale for large values of (f>. As 
the field turns round, it can drive a phase of inflation and, 
if the potential exhibits a sufficiently steep middle section 
around (f> ~ 0, reheating may proceed through gravitational 
particle production. Consequently, the field may survive until 
the present epoch, where it can act as the source of dark 
energy by slowly rolling along the tail towards 4> — * — °°- 



responds to the unstable saddle point that is also present 
in GR. The important characteristic of this solution is 
that any perturbations, no matter how small, necessar- 
ily force the universe to deviate exponentially from the 
static state. The second solution LS is always in the 
semi-classical domain and corresponds to a centre. This 
is a solution made possible by LQC effects. We have 
shown that it exists for a wide range of values of the po- 
tential, including positive, zero and negative values. Its 
importance lies in the fact that it allows a universe that 
is slightly displaced from the static state to oscillate in 
the neighbourhood of the static solution for an arbitrarily 
long time. 

We have exploited these characteristics to develop a 
working scenario of the emerging inflationary universe, 
in which a past-eternal, cyclic cosmology eventually en- 
ters a phase where the symmetry of the oscillations is 
broken by the scalar field potential, thereby leading in 
principle to a phase of successful slow-roll inflation. The 
mechanism that enables the universe to emerge depends 
very weakly on the form of the potential, and requires 
only that it asymptotes to a constant value at — > — oo 
and grows in magnitude at larger 4> in order to break 
the cycles. The asymptotic value of the potential can be 
either positive, zero or negative. 

The above emergent scenario has a further advantage 
that the initial state of the universe is set in the more 



natural semi-classical regime, rather than the classical 
arena of GR. Nevertheless, an important question that 
arises is the likelihood of these initial conditions within a 
more general framework. In particular, there is the issue 
of why the scalar field should initially be located in the 
asymptotic low-energy region of the potential. Although 
a detailed study of such questions is beyond the scope of 
the present paper, it has been argued |34| that for the 
case of a constant potential, the wavefunction in LQC 
most closely resembles the Hartle-Hawking no-boundary 
wavefunction j3^ . More specifically, the difference equa- 
tion in LQC requires the wavefunction to tend to zero 
near to the classical singularity |36| , an d in this sense re- 
sembles DeWitt's initial condition [33. However, within 
the context of solutions to the Wheeler-DeWitt equa- 
tion, requiring the wavefunction to be bounded as a — > 
selects the exponentially increasing WKB mode [34| and 
this corresponds to the no-boundary wavefunction. 

This is of interest since the square of the wavefunc- 
tion in quantum cosmology is interpreted as the proba- 
bility distribution for initial values of the scalar field in 
an ensemble of universes. For the no-boundary bound- 
ary condition, the probability, P oc exp(3/[8nlpV((j))]), 
is peaked at V{(f)) — 0, thereby suggesting that the no- 
boundary proposal disfavours inflation. In the present 
context, however, this implies that the most natural ini- 
tial condition for the field is to be located either at the 
minimum of the potential, or in the case where the poten- 
tial has no minimum, at (f) = —00. It would be interesting 
to explore this possibility further. In particular, previous 
analyses have so far neglected kinetic terms in the matter 
Hamiltonian and these may be important in the emerg- 
ing universe scenario. There is also the related question 
of whether the field moves from left to right initially. 

In general, we have found that the onset of slow-roll in- 
flation will occur at a relatively high energy scale, unless 
the quantization parameter j is extremely large. This in- 
dicates that a large amount of slow-roll inflation should 
arise, at least for a wide class of smoothly varying po- 
tentials, and it is expected, therefore, that the density of 
the present-day universe should be exponentially close 
to the critical density, Qo = 1 + e, where e <C 1. In prin- 
ciple, therefore, the emergent scenario we have proposed 
could be ruled out if a significant detection of spatial 
curvature is ultimately reported by future cosmological 
observations. 

It is also worth remarking that if the field is pushed 
sufficiently far up the potential, perhaps as high as the 
Planck scale, the inflating universe may enter a phase 
of eternal self-reproduction, where quantum fluctuations 
in the inflaton become more important than its classical 
dynamics 0. 

Finally, there is an interesting symmetry in the emerg- 
ing quintessential scenario between the initial and final 
states of the universe. Although the size of the universe 
differs by many orders of magnitude, the field evolves 
along the tail of the potential at both early and late times, 
4>(t — > —00) = <p{t — > +00), but with its kinetic energy 
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having changed sign. This implies that a reconstruction Acknowledgments 
of the dark energy equation of state at the present epoch 
could yield direct observational insight into the nature of 

the pre-inflationary potential in this scenario. DJM is supported by PPARC. 
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